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ABSTRACT 

We examine the fraction of massive asymptotic giant branch (AGB) stars remaining bound in their parent 
star clusters and the effect of irradiation of these stars by intracluster ultraviolet (UV) field. We employ a set of 
N-body models of dynamical evolution of star clusters rotating in a galactic potential at the solar galactocentric 
radius. The cluster models are combined with stellar evolution formulae, a library of stellar spectra, and 
simple models for SiO photodissociation in circumstellar environment (CSE). The initial stellar masses of 
clusters are varied from 50M 0 to 1O 5 M 0 . Results derived for individual clusters are combined using a mass 
distribution function for young star clusters. We find that about 30% of massive AGB stars initially born in 
clusters become members of the field population, while the rest evolves in star clusters. They are irradiated by 
strong intracluster UV radiation resulting in the decrease of the photodissociation radius of SiO molecules, in 
many stars down to the dust formation zone. In absence of dust shielding, the UV photons penetrate in the CSE 
deeper than 10R* in 64% and deeper than 2R, in 42% of all massive AGB stars. If this suppresses following 
dust formation, the current injection rate of silicate dust from AGB stars in the local Galaxy decreases from 
2.2 x 10 4 Mr 0 pc -2 Gyr 1 to 1.8 x 10 4 M ro pW 2 Gyr 1 at most. A lower revised value of 40% for the expected 
fraction of presolar silicate grains from massive AGB stars is still high to explain the non-detection of these 
grains in meteorites. 

Subject headings: stars: AGB and post-AGB stars: winds, outflows galaxies: star clusters: general (Galaxy:) 
solar neighborhood meteorites, meteors, meteoroids 


1. INTRODUCTION 

Circumstellar shells of low- and intermediate-mass stars 
(0.8 M q < M < 8M 0 4 ) at the thermally-pulsing AGB are 
known sites of dust condensation. Some grains from AGB 
stars, which finished their interstellar journey in material 
forming the solar system 4.6 Gyr ago, have been preserved 
in meteorites. Presolar origin of these grains (initial masses 
and metallicities of their parent stars) can be ascertained from 
anomalous isotopic ratios of the major and many trace ele¬ 
ments characterising stellar nucleosynthesis (Hoppe & Zinner 
2000; Nittler 2003; Dorschner 2010; Zinner 2014). 

Silicates and oxide grains are condensed in oxygen-rich 
stellar winds of AGB stars. The third dredge-up process mixes 
the carbon-rich ashes from He nuclear burning to the convec¬ 
tive envelope turning a star into a carbon star. This process 
becomes more efficient with increasing stellar mass, there¬ 
fore only low-mass AGB stars in the mass range of 1-1.5 M G 
are prominent sources of oxygen-rich dust. However, massive 
AGB stars with initial mass from 3.5-5 M 0 up to 8 M 0 , de¬ 
pending on model details and metallicity (Busso et al. 1999; 
Marigo 2007), convert the dredged-up carbon into I4 N as a re¬ 
sult of the hot-bottom burning process. There is recent spec¬ 
troscopic evidence that hot-bottom burning is active from the 
first thermal pulses through their final superwind phase (e.g., 
Garcia-Hernandez et al. 2013; Justtanont et al. 2012). 

Theoretical models of dust condensation in stellar winds of 


1 Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, D- 
86748 Garching, Germany 

2 Max Planck Institute for Astronomy, Konigstuhl 17, D-69117 Heidel¬ 
berg, Germany 

3 Institute for Astrophysics, University of Vienna, Tiirkenschanzstrasse 
17. A-l 180 Vienna 

4 We do not consider super AGB stars and assume that the upper limit is de¬ 
termined by the maximum mass of stars that develop an electron-degenerate 
CO core and end their life as white dwarfs (but see Doherty et al. 2013). 


evolved stars indicate that massive AGB stars should produce 
significant amount of oxygen-rich dust (e.g. Gail 2010; Fer- 
rarotti & Gail 2006; Ventura et al. 2012; Nanni et al. 2013). 
Models of the lifecycle of grains of different origins show that 
the oxygen-rich grains from massive AGB stars are expected 
to be an abundant component of the presolar grain popula¬ 
tion comparable to that from low-mass stars (Gail et al. 2009). 
The fact that grains are efficiently condensed in shells of AGB 
stars experiencing hot-bottom burning is supported by spec¬ 
troscopic observations of oxygen-rich dust-enshrouded stars 
(e.g.. Van Loon et al. 2005a). The chemical composition of 
these stars is altered by the hot-bottom burning process lead¬ 
ing to subsolar values of 12 C/ I3 C and l8 0/ 17 0 isotope ratios 
(Boothroyd et al. 1995), which should be possible to detect in 
presolar grains. However, no oxygen-rich grain with signa¬ 
tures of hot-bottom burning has been conclusively identified 
in meteorites so far. The entire population of presolar oxygen- 
rich grains appears to originate from low-mass AGB stars 
(Vollmer et al. 2008; Nittler et al. 2008; Iliadis et al. 2008; Nit¬ 
tler 2009; Hoppe 2010; Palmerini et al. 2011). To shed light 
at the origin of this discrepancy between theoretical models 
and findings from presolar grain studies, we consider addi¬ 
tional processes that can potentially reduce the global dust 
input from intermediate-mass stars. 

The total contribution of massive AGB stars to the Galactic 
dust budget depends on the adopted dust yields, i.e. on the 
amount of dust condensed in the stellar wind over the whole 
AGB evolution of a star. Dust yields as a function of the initial 
stellar mass and metallicity have been derived over the past 
decade from extensive calculations combining stellar evolu¬ 
tion at AGB stars, stellar wind and dust condensation models 
(e.g., Ferrarotti & Gail 2006; Ventura et al. 2012; Nanni et al. 
2013). Such models consider single isolated AGB stars un¬ 
der standard interstellar conditions. Nucleation and growth 
of dust particles in outflows of these stars is assumed to be 
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independent of their environment. This is justified for the ir¬ 
radiation of CSE by the standard interstellar radiation field, 
because the dust formation zone is located much deeper in 
the shell than the photodissociation region (e.g.. Jura & Mor¬ 
ris 1981;Glassgoldetal. 1987; Mamon et al. 1987; Glassgold 
1996). This may not be the case for evolved stars in young star 
clusters. Unlike low-mass stars, intermediate-mass stars are 
more likely to be members of their parent star clusters during 
their AGB evolution. 

The majority of stars in the mass range 4M 0 < M < 8 M a is 
formed in massive star clusters, which survive cloud dispersal 
and become open clusters, as will be shown below. The disso¬ 
lution time of open star clusters exceeds ~ 200 Myr, the life¬ 
time of a 4 M 0 star (Kharchenko et al. 2009; Just et al. 2009; 
Baumgardt & Makino 2003). The circumstellar environment 
(CSE) of massive AGB stars in young star clusters are subject 
to irradiation by the intracluster UV field from main sequence 
cluster stars of spectral types B and early A. Strong UV irradi¬ 
ation of CSE may affect the dust condensation process. Beck 
et al. (1992) investigated dust condensation in stellar winds 
of stars with ionising radiation of chromospheric origin and 
found that UV photons are able to reduce and, for certain field 
strengths, completely suppress the nucleation process. If the 
amount of dust condensed in outflows of massive AGB stars 
in star clusters is affected by external UV radiation, the net 
dust input from these stars will be reduced. 

In order to assess the importance of cluster environment on 
CSE of massive AGB stars, first of all we need to estimate the 
fraction of intermediate stars in the local Galaxy that evolve 
in their parent star clusters. To this end, we employ numeri¬ 
cal simulations of the dynamical evolution of star clusters in 
a mass range of 5O-1O S M 0 coupled with stellar evolution pre¬ 
scriptions (Kharchenko et al. 2009; Just et al. 2009; Ernst et al. 
2010). Next, for the massive AGB stars that remain bound 
in the cluster, we calculate the ultraviolet radiation field cre¬ 
ated by the neighbouring main sequence cluster members as a 
function of cluster mass. For a simple estimate of the impact 
of the intracluster radiation field on the circumstellar shells of 
AGB stars in cluster, we consider photodissociation of the SiO 
molecule. The SiO molecules are the basic building block for 
silicates, the most widespread dust species in oxygen-rich en¬ 
vironment (Henning 2010; Gail & Sedlmayr 2014). They are 
formed in the stellar photosphere, and in the shell their abun¬ 
dance decreases due to two processes: depletion on grains in 
the dust condensation zone and then photodissociation caused 
by interstellar UV photons in the outer regions. We calculate 
the photodissociation radius of SiO molecule for each mas¬ 
sive AGB stars in the model clusters and compare these radii 
to the position of the dust formation zone. 

The structure of our paper is as follows. Section 2 describes 
a suite of dynamical models of star cluster evolution, size- 
frequency distribution of star clusters and our main assump¬ 
tions in calculations of the UV fluxes irradiating circumstellar 
shells of massive AGB stars. Results of model calculations 
of dynamical evolution of model star clusters are presented in 
Sect. 3. We show that about 70 % of massive AGB stars ini¬ 
tially born in clusters remain cluster members. For these stars, 
we investigate the location of the photodissociation radii of 
SiO molecule with respect to the dust formation zone in their 
CSE. Dependence of our main results on the initial mass- 
radius relation for star clusters is presented in Sect. 4. Dis¬ 
cussion and concluding remarks are given in Sect. 5. 

2. MODEL 


2.1. Simulations of star cluster evolution 

The dynamical evolution of star clusters (SCs) is modelled 
as pure N-body systems of gravitationally bound stars rotating 
in a galactic potential at the solar galactocentric radius of 8.5 
kpc. The simulations were carried out with the direct /V-body 
code 0GRAPE (Harfst et al. 2007) on high-performance com¬ 
puting clusters with accelerator hardware GRAPE or GPU 
selected by appropriate libraries. These models were intro¬ 
duced and applied to individual SCs in a series of papers 
(Kharchenko et al. 2009; Just et al. 2009; Ernst et al. 2010). In 
the following, we briefly summarise the main model parame¬ 
ters used in our calculations and refer to the original publica¬ 
tions for a detailed description. 

2.1.1. Initial mass-radius relation 

We consider the evolution of low- and intermediate-mass 
stars in clusters, after the last massive stars ended their life 
as supernova, corresponding to the cluster ages of 40 Myr. 
Initially star clusters are embedded in their parent molecular 
clouds which are completely cleared off by stellar feedback by 
the onset of our simulations. Modelling of embedded cluster 
evolution and their emergence from parent molecular clouds 
requires including of gas dynamics or special algorithms ac¬ 
counting for the influence of gas on star dynamics as dis¬ 
cussed in Pelupessy & Portegies Zwart (2012) and Baumgardt 
& Kroupa (2007). These authors showed that the dynamical 
friction drives early mass segregation resulting in preferential 
retention of stars with mass above 2M 0 in clusters. Therefore 
we neglect possible ejection of stars during the early cluster 
evolution before the onset of our simulations. 

For the initial distribution of stars in the clusters, we adopt 
commonly used not-rotated King models (King 1966) with a 
central potential of Wo = 6.0 (Einsel & Spurzem 1999). In 
this case, the initial mass distribution is given by three pa¬ 
rameters: the initial cluster mass, concentration parameter 
and scale radius. Observed shapes of SCs as well as their 
sizes for a given cluster mass vary in a large range and can 
poorly constrain the initial conditions for simulations (Larsen 
2004; Portegies Zwart et al. 2010). This is due to the fact that 
they are influenced by a number of physical processes such as 
cluster formation, internal processes (self-gravitation, binary 
fraction, rotation), external forces (tidal field of the galaxy), 
encounters with giant molecular clouds (Scheepmaker et al. 
2007; Piskunov et al. 2008b; Kharchenko et al. 2009). Im¬ 
portance of these processes for the initial conditions of SC 
models is reviewed by Kroupa (2008). To minimize the num¬ 
ber of free parameters, we adopt a relation between the initial 
mass and radius 

R c i » 50 sJm c i/10 6 M q (pc) (1) 

derived from observations of molecular clouds and clumps in 
the Milky Way (Rivolo & Solomon 1988; Larson 1981). This 
mass-radius relation (MRR) has been used to set up the initial 
conditions for simulations of open clusters in the Milky Way 
(Adams et al. 2006; Fatuzzo & Adams 2008; Kharchenko 
et al. 2009). It has been also applied in theoretical studies 
of tidal tail clumps of SCs (Just et al. 2009) and for calibra¬ 
tion of radii and masses of observed clusters with simulations 
(Ernst et al. 2010). The initial half-mass radius of each cluster 
was adjusted so that the half-mass radius of a King model was 
initially equal to 80 per cent of a scaling radius and inside of 
which 60 per cent of cluster mass is enclosed. 

Given that there is no strong observational constraints or a 
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TABLE 1 

Initial parameters of the N-body models 


N, 

ALm 

M d 

(M 0 ) 

R hm (MRR) 
(pc) 

Rhm (TRR) 
(pc) 

^4</m/M q <8 

(V2«m/M 0 <8 

109 

4 

50 

0.25 

0.75 

1,2,1,1 

3,3,4,5 

153 

4 

70 

0.35 

0.84 

1,2,2,2 

5,5,6,6 

218 

4 

10 2 

0.44 

0.95 

1,2,2,2 

7,4,7,8 

437 

4 

2 x 10 2 

0.54 

1.20 

3,2,5,5 

15,9,17,15 

1094 

4 

5 x 10 2 

0.85 

1.63 

11,12,10,9 

30,28,40,27 

1532 

4 

7 x 10 2 

1.04 

1.82 

16,21,13,12 

50,51,47,54 

2189 

4 

10 3 

1.29 

2.05 

19,25,21,23 

60,69,66,77 

4379 

4 

2 x 10 3 

1.81 

2.58 

36,49,48,46 

130,155,140,142 

10947 

4 

5 x 10 3 

2.85 

3.50 

100,118,105,111 

364,373,345.351 

15326 

4 

7 x 10 3 

3.30 

3.92 

164,132,142,149 

519.510,492.492 

21895 

1 

10 4 

3.99 

4.41 

203 

716 

43791 

1 

2 x 10 4 

5.69 

5.56 

415 

1454 

109476 

1 

5 x 10 4 

8.96 

7.55 

1025 

3620 

153268 

1 

7 x 10 4 

10.62 

8.44 

1482 

5096 

218955 

1 

10 5 

12.66 

9.51 

2059 

7272 


Notes. The first column gives the total number of stars, followed by the number of simulation runs, the initial cluster mass, the initial half-mass radius for the 
mass-radius relation and that for the tidal radius-mass relation. The last two columns give the number of massive AGB stars and the number of irradiating stars 
in each simulation. 


universal relation for the initial conditions of simulated clus¬ 
ters, we test the dependence of our results on the adopted ini¬ 
tial cluster mass-radius relation and perform two additional 
sets of simulations. We follow an alternative method to set up 
the initial conditions assuming that the initial cluster radius is 
determined by the tidal field of the Galaxy and is given by (In- 
nanen et al. 1983; Baumgardt & Makino 2003; Lamers et al. 
2013) 

fltidai = 0 0.5GM cl /V*) 1/3 R 2 c f, (2) 

where M c \ is the cluster mass, Rq the galactocentric distance 
of the cluster and Vq the circular velocity of the Galaxy. This 
set of models is denoted as TRR. The half-mass radii for TRR 
and MRR models are shown in Table 1. 

Additionally, we perform a simulation run with the constant 
initial radius for all clusters (CRR models) to include the case 
of a weak mass-radius relation suggested by some observa¬ 
tions of galactic and extragalactic star clusters (Larsen 2004; 
Scheepmaker et al. 2007). We adopt a value of the initial half¬ 
mass radius of 2 pc (Lamers et al. 2013), consistent with the 
range of values from observations. In the following, we use 
the MRR models as the reference models and refer to them in 
the analysis of the results, unless stated otherwise. Compari¬ 
son of the main results for the TRR and CRR models with the 
reference models are given in Sect. 4. 

2.1.2. Stellar mass distribution in SCs 

The stellar mass in clusters is distributed following the ini¬ 
tial mass function (IMF) from Kroupa (2002) where only stars 
in the mass range of 0.08 < in <8.0 M 0 are studied. Stel¬ 
lar evolution of stars in clusters (stellar luminosities, radii, 
effective temperatures at various evolutionary phases) is de¬ 
scribed accordingly to the analytical formulas from Hurley 
et al. (2000). Binaries as well as multiple star systems were 
not included in our study. Possible encounters of SCs with gi¬ 
ant molecular clouds as well as passages of spiral arm or disk 
shocking were neglected. 

The initial mutual distribution of stars of different mass in 
clusters is still debated. There are many observational stud¬ 
ies in support of the initial mass segregation (IMS) of stellar 


mass, i.e. concentration of more massive stars in the inner 
cluster region (see Lada & Lada 2003; Portegies Zwart et al. 
2010, and references therein). However, observational evi¬ 
dence of mass segregation was questioned in Ascenso et al. 
(2009), who pointed at difficulties to differentiate between 
segregation and sample incompleteness effects. There are ob¬ 
servations that indicate that very low-mass stars associated 
with a cluster are distributed homogeneously in a volume that 
is much larger than the core of a cluster (Kumar & Schmeja 
2007). Given these uncertainties in the initial mass distribu¬ 
tion in star clusters, we run two sets of models: with and 
without initial mass segregation to study how it affects the 
evolution of massive AGB stars. Initial conditions for clusters 
with segregation are derived accordingly to the procedure de¬ 
scribed by Baumgardt et al. (2008). It invokes a random num¬ 
ber generator to draw distributions of initial stellar masses, 
positions, and velocities of stars in clusters. 

2.1.3. Initial parameters of model clusters 

We consider the dynamical evolution of 15 clusters with 
and without initial mass segregation resulting in 30 setups of 
simulation models. Initial masses, radii and numbers of stars 
for each cluster are listed in Table 1 . All clusters are moving 
on circular orbits in the Galaxy represented by an analytic 
background potential described by Plummer-Kuzmin mod¬ 
els (Miyamoto & Nagai 1975) with circular velocity Vq = 
233 km s _1 . Model SCs are assumed to have solar metallicity. 
For small clusters, the number of stars within a given mass 
range depends on randomized discretization of the initial stel¬ 
lar mass over individual stars. Therefore, we generate sev¬ 
eral sets of initial conditions for star clusters with the initial 
masses M c i < 7 x 10 3 M G resulting in 90 simulation models in 
total for one set. The number of simulation runs for each clus¬ 
ter model is given in the second column of Table 1. The table 
also shows the total numbers of stars with mass m > 4M G , 
which experience hot bottom burning during TP-AGB evolu¬ 
tion, and with m > 2M G , which dominate the intracluster UV 
radiation field, in each simulation run. Accordingly to Weid- 
ner et al. (2013) (see also Kroupa et al. 2013) the lowest mass 
of a star cluster to host at least one massive TP-AGB star is 
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~ 50 M q . Therefore this value is adopted here for M™ ln . We 
adopt an upper mass limit of model clusters M™ ax = 10 5 M o , 
appropriate for the Milky Way (Rruijssen 2014, and refer¬ 
ences therein). Star clusters of higher mass in the Galaxy 
belong to metal-poor >10 Gyr old globular cluster systems. 
Given short lifetimes of grains in the ISM (Slavin et al. 2015), 
dust from massive AGB stars from globular clusters could not 
survive until the solar system formation and is excluded from 
our consideration. 

Relevant information about all stars from the simulation 
runs (e.g., their coordinates, masses, luminosities, tempera¬ 
tures) is stored in snapshots with a time step of ~ 0.1 Myr. 
It is thus ensured that each massive TP-AGB star appears in 
several snapshots. 

2.2. Fraction of massive AGB stars in clusters 

In order to combine the results derived for individual star 
clusters of different masses to estimate the total fraction of 
massive AGB stars evolving in clusters, we need to know 
the cluster mass distribution function for the onset of simula¬ 
tions corresponding to the cluster age of 40 Myr. We assume 
that most of stars are formed in clusters embedded within gi¬ 
ant molecular clouds (Lada & Lada 2003; but see also Rrui¬ 
jssen 2012; Adamo et al. 2015). The distribution function 
of embedded star clusters can be measured observationally 
by infrared star counts. However, only the most massive 
(M c i > 5OOM 0 ) embedded clusters survive and become sta¬ 
ble, open clusters (Lada & Lada 2003; Bonatto & Bica 201 1). 
Since by the age of 40 Myr all clusters should emerge from 
their birth places, we can adopt the mass distribution func¬ 
tion of embedded star clusters, with a correction for the infant 
mortality of low-mass star clusters as described below. 

2.2.1. Embedded star cluster mass function 

Observations indicate that the mass-frequency distribution 
of young embedded star clusters dn em \,(M c \)/dM c \ follows 
a universal power law (Kennicutt et al. 1989; McKee & 
Williams 1997; Lada & Lada 2003; Larsen 2010; Porte- 
gies Zwart et al. 2010): 



log(M d , M 0 ) 

Fig. 1.— Cumulative distribution of progenitors of massive TP-AGB stars 
(in all clusters with mass below M c \) normalised to the total number of these 
stars in all clusters calculated using the embedded cluster mass function. The 
shaded area marks the initial mass of clusters, which emerge mostly unbound 
from their parental molecular clouds, accordingly to Lada & Lada (2003). 

Figure 1 shows the cumulative distribution of stars with ini¬ 
tial masses 4M 0 < M < 8M G , which are progenitors of mas¬ 
sive TP-AGB stars, as a function of initial cluster mass for the 
distribution function of young clusters given by Eq. (3). It is 
derived for the adopted Kroupa IMF and cluster mass range. 

Because of gradual dissolution of SCs, the actual values 
of bound AGB stars A^ gb (M c i) derived from simulations are 
lower. Additionally, the mass distribution function of embed¬ 
ded star clusters has to be corrected for the dissolution of low- 
mass clusters upon dispersal of molecular clouds as described 
below. 

2.2.2. Dissolution of low-mass clusters 


dn emb (M cl )/dM cl oc M c f, M” m < M cI < (3) 

where dn e mbCV c i) is the number of embedded star clusters 
with the masses between M c \ and M c \ + dM Q \, M™ n and M™ ax 
are the lowest and highest star cluster masses. As discussed 
above, we assume M™ n = 50M o and M™ ax = 10 5 M G , respec¬ 
tively. 

The total number of AGB stars formed in all clusters with 
the masses below M c \ is 

A^ B (M d ) = f d N AG B(M) dnc ^ M) dM. (4) 

in ClM 

cl 

where Nagb(M) is the initial number of AGB stars in a cluster 
with the mass M. 

The cumulative number of AGB stars formed in a cluster 
with mass M < M c \ normalised to the total number of these 
stars is easily derived from Eq. (4) 


r* A 

JM r 

f 


«<*«= 4 A, "“ (M)d w 5dM / 

N ACB m dnemb . (M) dM. 


(5) 


dM 


Numerical models in this work consider evolution of clus¬ 
ters shortly after emergence from molecular clouds, there¬ 
fore we include two sources of the field population of mas¬ 
sive AGB stars separately: (i) members of low-mass clusters 
(M c i < 500M G ) which did not survive dispersal of molecular 
clouds, and (ii) stars lost by bound clusters in the process of 
their dynamical evolution. The latter is derived from our dy¬ 
namical N-body models as a function of initial cluster mass. 

The fraction of star clusters that survive as bound systems 
up to Pleiades age was assessed by Lada & Lada (2003) by 
comparison of the numbers of embedded and open star clus¬ 
ters from observations. They found that majority (90-95%) 
of embedded star clusters with initial mass M c \ < 500M 0 
emerge from molecular clouds as unbound systems and con¬ 
tribute their members to the field population. The contribution 
of clusters in this mass range to the total number of progen¬ 
itors of massive TP-AGB stars is about 30%, while about a 
half of progenitors of massive AGB stars are formed in clus¬ 
ters with the mass above 2 x 10 3 M G (Fig. 1). 

Dynamical models which consider embedded evolution of 
star clusters in molecular clouds find that the survival of em¬ 
bedded star clusters depends on many physical parameters 
such as gas expulsion timescale, star formation efficiency, and 
impact of an external tidal field (Parmentier & Baumgardt 
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2012; Parmentier 2011). Nevertheless, the slopes of the mass 
functions of embedded or young star cluster and gas-free clus¬ 
ters are identical (Dowell et al. 2008; Oey et al. 2004; Lada 
& Lada 2003). For a simple estimate of the contribution of 
clusters M Q \ < 500M 0 to the field population, we assume that 
initially these SCs were formed accordingly to the mass distri¬ 
bution function dn em \,(M)/dM given by Eq (3) as their higher 
mass counterparts, but 90% of these clusters do not survive 
cloud dispersal based on the lower limit derived by Lada & 
Lada (2003) and more recently by Bonatto & Bica (2011). 
Therefore, for the onset of our simulation, we adopt the mass 
distribution function of young star clusters dzz colT (M)/dM cor¬ 
rected for the dissolution of low-mass clusters by multiplica¬ 
tion of dn em b(M)/dM by 0.1 for M c \ < 500 M Q . 

The total number of bound AGB stars is 


N 1 


-bound 

AGB 


r M ™ 

Jmt 




d^corr( T/) 

d M 


d M 


( 6 ) 


where N b v]ii (M) is the number of bound stars in a SC with 
mass M derived from numerical simulations of dynamical SC 
evolution. The total number of bound AGB stars residing in 
star clusters with M < M c i relatively to the number of all AGB 
stars therefore equates to 


/AGB(^d) = 


yv tot 

JV AGB 


nM c \ 

I J 

JmT 1 




d^corrl '1'/ ) 

d M 


d M, (7) 


where N^, B is given by Eq. (4) for M c \ = M™ x 


2.3. Intracluster UV radiation field 

Circumstellar envelopes of evolved stars in young star clus¬ 
ters are irradiated by UV radiation emitted by the main se¬ 
quence cluster members of spectral type A and late B. In 
older clusters, thermal radiation from white dwarfs becomes 
an important source of UV photons, which is probably respon¬ 
sible for the ionisation and removal of intracluster medium 
(McDonald & Zijlstra 2015). The strength of the intracluster 
UV radiation field due to main sequence stars is expected to 
quickly decrease over a few hundred Myr as the cluster turn 
off point shifts towards stars of lower mass, since amount of 
emitted UV radiation is a strong function of stellar mass. It 
is illustrated in Fig. 2 showing the total number of ionising 
photons for H I and C I atoms emitted per unit time Qq as 
a function of initial mass for intermediate-mass stars. Qq is 
derived by integration of the stellar surface flux F v 


Qo 



4 nR;F v 
hv 


dv. 


( 8 ) 


where IF is the stellar radius, and v max and v m ; n correspond 
to the wavelength of 91.178 nm and 240 nm (5.17 eV) for 
H I and to 91.178 nm and 110.11 nm (11.26 eV) for C I , 
respectively. 

Intermediate-mass main sequence stars in clusters are capa¬ 
ble to ionise the intracluster medium. Their radiation is hardly 
attenuated within the cluster owing to the observed lack of 
gas and dust in the intracluster medium (Bastian & Strader 
2014). Figure 2 shows the radii of Stromgren spheres, i.e. 
zones of ionised H II and C II around the central star, for the 
homogeneous intracluster gas density of 1 cm~ 3 , similar to 
the density of the ambient ISM. Although the sizes of ionised 
regions around intermediate-mass stars are small compared to 
those around massive stars, they are comparable to the sizes 
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Fig. 2.— Number of ionising photons for H and C atoms emitted per unit 
time as a function of initial stellar mass (left axis) and corresponding radii of 
ionised zones (right axis) calculated for the value of ambient gas density of 
n\\ = 1 cm -3 . The second axis on the right shows log R s in pc. 


of low-mass clusters (Table 1). After the onset of mass loss 
from evolved stars in clusters, the density distribution in the 
vicinity of these stars will follow F 1 radial dependence and 
only outer layers will be ionised by the extreme UV radiation. 

The strength of the UV field is customary characterised rel¬ 
atively to the standard ISRF (Rollig et al. 2007) 

/'»240nm r»240nm 

U = Au. u dAl Auj SRF dA, (9) 

v/91.2nm J91.2nm 

where uj = An/cJ,\ the spectral photon energy density, m*j SRF 
is the the spectral photon energy density of the interstellar ra¬ 
diation field (Draine 1978). The mean radiation field intensity 
in the location of an AGB star is approximated as 

«MS «MS ft! 

4 ttJ.i = ^ I\ t + AnJ isrf = ^ + AnJ isrf (10) 

;= l i=i r i 


where zzms is the current number of main sequence stars with 
m, > 2M 0 in a cluster, I\, is the incident intensity from zth 
star, Rj is its photospheric radius, and r, is the distance from 
the central AGB star to zth star, and Q' is the solid angle 
subtended from the zth star. The last term on the right side 
of equation is the contribution from the interstellar radiation 
field. t is derived from the Eddington fluxes H v provided 
by a stellar spectral library 


h = 


Fa 

n 



(ID 


The stellar fluxes are assigned to cluster stars accordingly to 
their temperature and surface gravity stored in the snapshots. 
We employ the BaSeL 3.01 semi-empirical library of stellar 
spectra for wavelength from 9 nm to 160 ftm (Lejeune et al. 
1997, 1998; Westera et al. 2002). 


2.4. Photodissociation radius ofSiO molecules 

SiO molecule is considered a prerequisite for formation of 
quartz and various silicate types of grains (Henning 2010; 
Gail & Sedlmayr 2014). It is well established by photochemi¬ 
cal models and supported by observations that the abundance 
of SiO molecules decreases in the CSE of AGB stars with 
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radius by two processes: depletion on grains in dust con¬ 
densation zone and dissociation by interstellar UV photons 
in outer regions (Jura & Morris 1981; Huggins & Glassgold 
1982a,b; Gonzalez Delgado et al. 2003). Such photochemical 
models usually assume a homogeneous dust distribution with 
dust attenuation properties similar to those in the interstellar 
medium. In this case, dust shielding plays a crucial role in 
the chemistry of CSE, being able to protect molecules from 
the interstellar UV photons up to distances of 10 16 - 10 17 cm 
(Glassgold 1996). In this work, we evaluate how deep the UV 
photons can penetrate in the expanding shell of an AGB star, 
which is irradiated by the enhanced UV radiation prior to the 
onset of efficient dust condensation. We consider the sizes of 
the SiO envelopes around AGB stars in the model clusters and 
compare them to those irradiated by the ambient ISRF, and to 
the position of the dust formation zone. 

For a rough estimate, we adopt a toy model of the photodis¬ 
sociation of SiO molecule in a stationary spherically symmet¬ 
ric outflow, where the gas has been accelerated to terminal 
velocity v e . In absence of other sinks and sources, the SiO 
fractional abundance /s;o relatively to H 2 changes due to pho¬ 
todissociation by UV photons as (Jura & Morris 1981; Hug¬ 
gins & Glassgold 1982b) 


d/sio 

dr 


t pd 

^SiO 


■ exp 



( 12 ) 


where dsio is the shielding distance corresponding to the op¬ 
tical depth of 1 and k^. () is the unattenuated photodissociation 
rate of SiO molecules. For a monosize grain distribution, the 
shielding distance is 


dsiO - 1-4 


3(6/fl)sio 

4psio 


M d 

X 4nv e ’ 


(13) 


where Q is the dust absorption efficiency, a is the grain size, 
Psio is the density of the solid material, M d is the dust mass- 
loss rate. The grain drift relatively to the gas was neglected. 
The value of dust mass loss is M d = DM, where D is the dust- 
to-gas ratio in the stellar wind and M is the gas mass-loss rate 
from Eq. (16). 

The size of the SiO envelopes derived from observations 
agrees well with the photodissociation radius defined as the 
distance from the central star where the SiO abundance is de¬ 
creased by e times (Gonzalez Delgado et al. 2003). The equa¬ 
tion for the photodissociation radius is derived from Eq. (12) 
(Jura & Morris 1981) 


' P ~ n d ' 

klt 0 E 2 (d si0 /r p ) 

where Ei(x) denotes the exponential integral. Equation (14) 
is solved numerically for all AGB stars in the model clusters. 
In the dust-free case, the photodissociation radius is simply 
given by r p = v e lk pd Q . 

The unattenuated photodissociation rate of SiO molecule 
through line absorption from the lower level 1 into an upper 
state level u is (van Dishoeck et al. 2006) 


k 


pd 

SiO 


ne 2 4nJ Y (v uI ) 

m e c hvui 


05) 


where /„/ is the line absorption oscillator strength, //„ is the 
dissociation efficiency of the upper state u, which lies be¬ 
tween 0 and 1. Assumption of p,, - 1 is reasonable for SiO 


molecule. Mutual shielding and self-shielding are the factors 
that can reduce the photodissociation rate of such molecules 
as CO, H 2 . Mutual shielding is probably not important for 
SiO (E. van Dishoek, private communication). The total pho¬ 
todissociation rate is computed by summing over all lines. We 
include absorption to the 3 1 X' and 2, 3, 4, and 5 'n states for 
SiO molecule with the oscillator strength / = 0.10, 0.32, 0.03, 
0.11, and 0.10, respectively, taken from the photodissociation 
database 5 (van Dishoeck et al. 2006). The wavelengths for 
absorption lines are 1011, 1058, 1063, 1140, and 1378 A. 

The mean intensity J v is derived from Eq. (10), for the in¬ 
tracluster radiation field, and is taken from work by Mathis 
etal. (1983), for the ISRF. 


2.4.1. Wind model parameters 

We focus on dust formation by massive AGB stars with hot- 
bottom burning, therefore we consider only an oxygen-rich 
dust mixture in circumstellar envelopes. Although the dust 
composition in envelopes of M-stars depends on the initial 
stellar mass and mass-loss rates, for simplicity we assume a 
fixed silicate composition with density p d = 3.3 gem 3 . We 
adopt the following values of the dust-to-gas ratio D in the 
wind: D — 0 corresponding to a dust-free case, and D = 
0.001, a typical value in outflows after the condensation zone 
derived from theoretical models of massive AGB stars. Grains 
are assumed to have a single size of 50 nm. The absorption 
efficiency Q for a grain size of a — 50 nm and A ~ 100 nm 
equates to 1. The optical depth at around 100 nm, calculated 
for the adopted grain size, is two times lower compared to that 
calculated for a power law distribution of grain sizes from 
work by Kim et al. (1994), which is commonly assumed in 
spectral energy distribution fitting. 

The mass-loss rate is assumed to be isotropic. We adopt the 
mass-loss rate for oxygen-rich AGB stars derived empirically 
by Van Loon et al. (2005a), which is more suitable for massive 
AGB stars than Reimer’s law 


M = 1.38 x 10" 11 



( 


Eeff V 603 

3500 K/ 


(16) 


where L and T cn - are the luminosity and effective temperatures 
of TP-AGB stars. The expansion velocity v e is assumed to be 
10 kms -1 . For Si abundance, the value esi = 4 x 10 7 is 
adopted. It is assumed that initially all Si atoms are bound in 
SiO molecules. 


3. RESULTS 

In this section, we present results of model evolution of 
massive AGB stars in star clusters of various initial mass (Ta¬ 
ble 1) rotating in the Galaxy at the solar galactocentric radius 
from simulation time t - 40 Myr to 220 Myr corresponding 
to the lifetimes of these stars. 

3.1. The fraction of massive AGB stars in clusters 

The relative fractions of bound and unbound massive AGB 
stars for each model cluster are displayed in Fig. 6. The frac¬ 
tion is derived from the analysis of binding energies of all 
massive AGB stars in the cluster for the entire simulation run. 
Most of unbound massive AGB stars originate from low-mass 
star clusters. While low-mass clusters with initial mass seg¬ 
regation eject a slightly higher number of massive AGB stars 

s http: //home. strw. leidenuniv. nl/-ewine/photo/ 
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Fig. 3.— Intracluster UV radiation field in the locations of massive AGB 
stars in star clusters with initial mass of 10 2 and ! O' /V/. 0 relatively to the 
standard interstellar radiation field calculated for each simulation snapshot. 
Different colors show four simulation runs. Top and bottom panels show 
models without and with initial mass segregation, respectively. 


compared to the models without IMS, it is the opposite for 
massive clusters (M c \ > 2 x 1O 3 M 0 ) with IMS keeping al¬ 
most all their massive AGB stars. The fractions of massive 
AGB stars lost by massive clusters modelled without IMS are 
0 . 2 - 6 %. 

Figure 7 shows the cumulative distribution of massive AGB 
stars calculated with account of the mass distribution function 
of star clusters as described in Sect. 2.2. Most of unbound pro¬ 
genitors of massive AGB stars that eventually become mem¬ 
bers of the field population originate from star clusters with 
low-mass clusters, while more massive clusters preserve their 
massive AGB stars (see Fig. 6). The majority of low-mass 
clusters are dissolved in the process of emergence from their 
parent molecular clouds before the onset of our simulations. 
It is taken into account in the same way for both types of ini¬ 
tial conditions, therefore the cumulative distributions derived 
for models with and without initial mass segregation appear 
very similar in Fig. 7. Our main finding is that 70% of mas¬ 
sive AGB stars at the solar galactocentric radius reside in star 
clusters. 30% of all massive AGB stars were lost from their 
clusters prior to the onset of AGB evolution and evolve in 
the field. They originate mainly from star clusters with initial 
mass M c i < 10 3 M 0 . In the following we focus on massive 
AGB stars in cluster environment and its effects on their cir- 
cumstellar envelopes. 

3.2. Intracluster UV radiation field 


Fig. 4.— Intracluster UV radiation field at the locations of massive AGB 
stars in star clusters with initial mass of 10 4 M e without and with initial 
mass segregation (top and bottom, respectively) calculated for each simu¬ 
lation snapshot. The flux is shown relatively to the standard interstellar ra¬ 
diation field. Red squares, green stars, blue crosses and black plus symbols 
show UV field in the location of AGB stars in the snapshots with one, two, 
three and four AGB stars, respectively. 

In the following, we investigate the strength of the intra¬ 
cluster radiation field, its time evolution and dependence on 
cluster mass using numerical simulations described in the 
Sect. 2.1. Figures 3-5 show the strength of unattenuated 
intracluster radiation field in the locations of massive AGB 
stars for star clusters with initial mass M c \ = 10 2 , 10 3 , 10 4 , 
and 10 5 M g . The intracluster radiation field is calculated for 
each snapshot by summing up the fluxes of photons in the en¬ 
ergy range of 5.6-13.598 eV from all cluster main sequence 
stars with mass m > 2M 0 . In order to compare the irradiation 
which an AGB star experiences in a star cluster to that in the 
ambient ISM, the fluxes are shown relatively to the standard 
interstellar radiation field (ISRF). We adopt the value of the 
photon flux /-j) raln e = 1.921 x 10 8 cm~ 2 s _1 (Draine 1978). 

Since the duration of thermally pulsing AGB stage is rel¬ 
atively short, low-mass clusters have no more than one star 
in this evolutionary phase at the same time (Fig. 3). A clus¬ 
ter with M c i = 10 4 M 0 has up to four TP-AGB stars at the 
same time, while a 10 5 M 0 cluster has generally more than 
four massive AGB stars. The UV fluxes at the location of 
multiple AGB stars are shown with different symbols depend¬ 
ing on their number in Figs. 4 and 5. 

Although there is a scatter in the values of the UV flux, 
figures 3 and 4 demonstrate that most of AGB stars are irra¬ 
diated by 10-100 times stronger UV field than the ISRF. The 
average irradiation flux decreases with time as stars of higher 
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Fig. 5.— The same as in Fig. 4 for M c \ = 10 5 M o . Grey circles indicate the 
cases, when there are more than four AGB stars in a snapshot. 




Fig. 7.— Cumulative number of massive AGB stars evolving bound in the 
clusters, unbound (ejected) and their sum (triangles, circles and squares, re¬ 
spectively) relatively to the total number of these stars as a function of initial 
cluster mass M c \. Solid and dashed lines indicate results for models with and 
without initial mass segregation, respectively. 


TABLE 2 

Number of massive AGB stars with R p d,siO < ^cond relatively to the 
NUMBER OF THESE STARS BOUND IN SCs 


IMS no IMS 


^cond/^* 

D = 0 

D = 10~ 4 

D = 0 

D = 10~ 4 

2. 

0.60 

0.04 

0.40 

0.05 

5. 

0.80 

0.25 

0.60 

0.14 

10. 

0.92 

0.64 

0.75 

0.41 


Notes. The first column gives the adopted values for the condensation radius 
/?cond relatively to the photospheric radius. It is followed by the relative 
numbers calculated for different values of the dust-to-gas ratio D in the CSE 
for models with and without IMS. 


TABLE 3 

Number of massive AGB stars with /? p d,siO < ^cond relatively to the 
NUMBER OF ALL MASSIVE AGB STARS 


IMS no IMS 


^cond /R* 

D = 0 

D = 10~ 4 

D = 0 

D = 10 

2. 

0.42 

0.03 

0.28 

0.03 

5. 

0.56 

0.18 

0.41 

0.09 

10. 

0.64 

0.45 

0.52 

0.29 


from their more extended stellar distribution. 


log(M d , M G ) 

Fig. 6.— Number of bound and unbound massive AGB stars (triangle and 
circle symbols, respectively) relatively to the number of these stars in each 
cluster as a function of initial cluster mass. Dashed and solid lines show star 
cluster models with and without IMS, respectively. 


mass go off the main sequence (Fig. 2). There is more than an 
order of magnitude spread in UV flux strength towards lower 
values for clusters without initial mass segregation resulting 


3.3. Photodissociation ofSiO molecules by intracluster field 

3.3.1. Photodissociation radii ofSiO molecules in cluster AGB 

stars 

We compute the photodissociation radii /? p d,sio of SiO en¬ 
velopes of bound AGB stars in all simulated clusters for two 
cases of UV irradiation: (i) by the intracluster radiation field 
in the location of each AGB star in model clusters, and (ii) by 
the standard isotropic ISRF, for comparison. Figure 8 shows 
the photodissociation radii for cluster masses of 10 3 , 10 4 , and 
10 5 M g for MRR models with and without IMS for dust-free 
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Fig. 8. — Photodissociation radius of SiO relative to photospheric radius R in dust-free shells of massive AGB stars in star clusters of 10 3 , 10 4 , and 10 5 (from 
top to bottom). Red symbols show results of calculations with dissociation only by the ISRF component, and blue symbols by the ISRF and intracluster UV 
photons. Left and right panels present results for N-body models with and without IMS, respectively. The grey area marks the dust formation zone. 


stellar wind. This case corresponds to the onset of dust con¬ 
densation process. Time variations of the photodissociation 
radii of AGB stars in the cluster environment is similar to 
those of the intracluster UV field strength (Fig. 3-5). The 
difference is that the former is determined by the intensity in 
SiO dissociation lines, while the latter characterises the con¬ 
tinuum UV radiation. 

During AGB phase, stars change their locations in the clus¬ 
ter, therefore they are irradiated by variable UV flux. From 
analysis of /? p d,sio in different snapshots for the same AGB 
stars, we find that the corresponding variations of /? p d,sio are 


typically less than 2 times. 

It is insightful to compare the derived photodissociation 
radii with the dust condensation zone in CSE. The location 
of this zone is a matter of debate. According to theoreti¬ 
cal models of dust condensation in outflows of M-stars, dust 
is formed at 5-10 photospheric stellar radii R t (e.g., Gail & 
Sedlmayr 1999; Ferrarotti & Gail 2003). Nanni et al. (2013) 
recently derived smaller radii for dust condensation of 2 - 3 R* 
using similar to the aforementioned models with an alterna¬ 
tive mechanism of dust destruction. A smaller condensation 
radius is also required by dynamical models of CSE of M-type 
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AGB stars (Hofner 2008). Observations of spatially resolved 
envelopes yield the values varying from a few to ~ 10 R* (e.g., 
Khourietal. 2014; Zhao-Geisleret al. 2012). We highlight the 
radii from 2 to 10 R t in Fig. 8 as the dust condensation zone. 

As expected, /s’pd.sio due to photodissociation by the inter¬ 
stellar UV photons lies outside of the dust formation zone in 
all clusters even without dust shielding. The cluster UV pho¬ 
tons penetrate much deeper in the CSE, reaching the dust con¬ 
densation regions in a significant fraction of stars. A common 
feature for all model clusters is that the initial segregation of 
stellar mass has strong impact on /x’ p d,sio owing to the higher 
stellar density of both irradiating and AGB stars in the central 
regions of clusters. 

Figure 9 illustrates the dependence of the SiO photodissoci¬ 
ation radii on the adopted dust-to-gas ratio D in the shell. The 
figure shows /fpd.sio for stars in the 5 X 10 4 M Q cluster and the 
dust-to-gas ratio D — 0, 10 4 , and 10 5 . Similarly to other 
models, the intracluster radiation field in this cluster has the 
largest impact on SiO in the shells of more massive AGB stars 
ending their life during the first 100 Myr. In the dust-free case, 
the 54% (18%) of stars have R p d.sio < 2 R t in models with¬ 
out) IMS. For D = 10 ~ 4 , /? p d,sio is noticeably pushed outside 
due to dust shielding. There are no stars with SiO envelope 
within 2 R*. Furthermore, 20% of stars have /? p d,sio <5 R* 
in models with IMS and none without IMS. The fraction in¬ 
creases to 60% ( 20 %) for our upper limit for the condensa¬ 
tion radius of 10 R* in the models with(-out) IMS. The ratio 
D = 10 ~ 3 provides sufficient attenuation to block dissociating 
photons from the dust formation zone. The intracluster radia¬ 
tion in this case does not affect the condensation process, but 
it reduces the extent of SiO envelopes from a few times up to 
an order of magnitude. This effect can be used in observations 
of massive AGB stars as an indicator of cluster membership. 

3.3.2. Fractions of stars affected by the intracluster field 

We calculate the total number of bound AGB stars with 
^pd.sio < Rcond for three values of the condensation radius 
f?cond = 2, 5, and 10 R, discussed above for each simulated 
cluster. These numbers are used to derive the total fractions 
of massive AGB stars affected by the intracluster UV field by 
integration over the entire cluster mass range with the cluster 
mass distribution function. Table 2 shows the resulting frac¬ 
tions relatively to the total number of bound massive AGB 
stars. It is derived using the total number of the bound stars 
N b ™* d given by Eq. ( 6 ) in the denominator of Eq. (7) instead 
of For N ^ GB , we take the number of bound stars satis¬ 

fying the condition A’| U i,sio < Rcond- The fractions of stars are 
shown for models with and without IMS for the dust-to-gas 
ratio values D — 0 and 10 4 . For the dust-free case, /f p d,sio 
in majority of AGB stars in clusters with IMS is deeper than 
the inner border of dust condensation zone. The fractions of 
AGB stars in clusters affected by intracluster field are lower 
in models without IMS, but still significant, 60% and 75% 
with 7? p d,sio <5 R, and Wpd.sio < 10 /?„, respectively. The 
numbers in Table 2 imply that even if minor refractory dust 
species are formed closer to the star resulting in the dust-to- 
gas ratio of < 10 ~ 4 , shielding provided by these species is not 
sufficient and the intracluster UV field nevertheless may affect 
dust condensation in a large fraction of stars. 

As discussed in Sect. 3.1, results of our dynamical models 
of star cluster evolution combined with the adopted distribu¬ 
tion of young star clusters imply that about 30% of all massive 
AGB stars evolve in isolation. To illustrate the relative impor¬ 


tance of cluster environment for the entire population of mas¬ 
sive AGB stars (bound and unbound), we show the numbers 
of AGB stars with /? pc i,sio < Rcond relatively to the number of 
all massive AGB stars in Table 3. According to our model 
calculations, in absence of dust shielding, the intracluster UV 
photons are able to penetrates in the CSE deeper than 10 R* 
in 52 - 64% and deeper than 2 R* in 28 - 42% of all massive 
AGB stars. The lower and higher values correspond to the 
models without and with IMS. 

3.4. Implications for the stardust input from AGB stars 

3.4.1. Revision of the model predictions for presolar grains 

In the following, we revise the model predictions for the rel¬ 
ative contribution of massive AGB stars to the presolar grain 
inventory carried out in our earlier work (Gail et al. 2009). 
Gail et al. (2009) put forward a model, which relates theoret¬ 
ical studies of dust condensation during the entire AGB stage 
and a model for chemical evolution of dust and gas in the solar 
neighbourhood to the contribution of AGB stellar populations 
as a function of their initial mass and metallicity. This model 
suggests that 40% and 60% of silicate grains in the stardust 
population at the instant of solar system formation originate 
from low- and high-mass AGB stars, respectively. 

In the present study, we can evaluate the maximum possible 
effect of cluster environment on the contribution from AGB 
stars by assuming that dust condensation is completely sup¬ 
pressed in massive AGB stars with R p d.sio < Rcond and adopt¬ 
ing the largest dust condensation radius R coa d = 1 OR , . We use 
the same dust yields as in our earlier works (e.g., Zhukovska 
et al. 2008; Gail et al. 2009; Zhukovska 2014) based on dust 
condensation models in O- and C-rich stellar winds of AGB 
stars. These yields are derived for single stars. There are 
observations that some binaries develop a common envelop 
which allows dust condensation process. Evolution of inter¬ 
acting binaries is poorly understood. There are presently no 
yields available for dust condensed in the common envelopes 
for the solar metallicity. Very recently, Zhu et al. (2015) es¬ 
timated the contribution of dust from binaries for the Large 
Magellanic Clouds and found that it is 4 times lower com¬ 
pared to the dust input from single stars using the dust yields 
from Zhukovska et al. (2008). 

We multiply the dust masses ejected by massive AGB stars 
by the factor 1 - 77 , where 77 is the fraction of stars with pho¬ 
todissociation radii R p d.sio < Rcond listed in Table 3. For 
I] = 0.64 taken for the IMS models, the corresponding re¬ 
vised value for the relative contribution of massive AGB stars 
to the O-rich presolar grains of the AGB origin equates to 0.4. 
We conclude that although including star cluster environment 
may reduce the relative contribution of massive AGB stars to 
the presolar grains, the lower limit for the mass fraction of 
40% is too high to help mitigate non-detection of grains from 
these stars in meteorites. 

3.4.2. Dust production rates 

AGB stars are the main stellar source of dust produc¬ 
tion in the Galaxy. The present day rate of dust injection 
from AGB stars at the solar galactocentric radius is 1.2 x 
10~ 3 M Q pc ~ 2 Gyr -1 , out of which 2.2 x 10~ 4 M G pc ~ 2 Gyr 1 is 
attributed to the silicate stardust injection rate, accordingly to 
the models published in Zhukovska et al. (2008). The con¬ 
tribution from massive AGB stars is 8 x 10~ 5 M o pc ~ 2 Gyr~' 
or 35% of silicate dust input from AGB stars. Assuming that 
dust formation is completely suppressed in the cluster stars 
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Fig. 9. — Photodissociation radius of SiO relative to photospheric radius R. in massive AGB stars in a star cluster with mass of 5 X 10 4 M Q and dust-to-gas ratios 
of 0, 10 4 , and 1 0 4 (top, middle and bottom panels, respectively). The meaning of the symbols is the same as in Fig. 8. Each second snapshot is shown to reduce 
crowding of the figure. 


with /? p d,sio < 10/?* as described above and multiplying the 
dust yields for m > 4M Q by the reduction factor 1 - j) of 0.46, 
we attain a lower limit for the current injection rate for AGB 
stars of 1.8 x 10~ 4 M G pc 2 Gyr’ 1 corrected for the fact that 
a substantial number of massive AGB stars remain bound in 
their parent star clusters. Thus including the effect of cluster 
environment in dust condensation models can reduce the to¬ 
tal injection rates of silicate dust from AGB stars by at most 
20 %. 


In any case, if the total silicate production rates from AGB 
stars are decreased owing to the effects mentioned above, it 
should not affect the amount or composition of the interstellar 
dust mixture, since the mass fraction of stardust from AGB 
stars in the Milky Way is only 2% and it is dominated by the 
carbonaceous grains. The majority of dust mass is formed by 
accretion in the ISM (e.g., Dwek 1998; Zhukovska et al. 2008; 
Draine 2009). 
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TABLE 4 

Number of massive AGB stars with Spd.siO < R C ond for the TRR models 



relatively to the number of these stars bound in SCs 

relatively to the number of all 

massive AGB stars 



IMS 


no IMS 


IMS 


no IMS 

^cond /R* 

D = 0 

D = 10- 4 

D = 0 

D = 10~ 4 

D = 0 

d = ur 4 

D = 0 

D = lO^ 4 

2. 

0.53 

0.02 

0.31 

0.03 

0.37 

0.02 

0.22 

0.02 

5. 

0.74 

0.21 

0.51 

0.09 

0.53 

0.15 

0.36 

0.06 

10. 

0.89 

0.57 

0.70 

0.34 

0.63 

0.40 

0.49 

0.24 





TABLE 5 







Number of massive AGB stars with R p d,siO 

< - ^cond 

FOR THE CRR 

MODELS 



relatively to the number of these stars bound in SCs 

relatively to the number of all 

massive AGB stars 



IMS 


no IMS 


IMS 


no IMS 

^cond/^* 

D = 0 

D = 10~ 4 

D = 0 

D = 10~ 4 

D = 0 

d = nr 4 

D = 0 

D = 10~ 4 

2. 

0.77 

0.05 

0.66 

0.02 

0.54 

0.03 

0.47 

0.02 

5. 

0.86 

0.60 

0.78 

0.42 

0.60 

0.42 

0.56 

0.30 

10. 

0.92 

0.79 

0.88 

0.69 

0.64 

0.55 

0.63 

0.49 


4. DEPENDENCE ON THE INITIAL CLUSTER 
MASS-RADIUS RELATION 

In this section, we discuss how alternative relations between 
the initial mass and radius of model clusters based on the TRR 
and CRR models affects our main conclusions. All initial pa¬ 
rameters for these models are the same as for the MRR set of 
models (Table 1) with exception of the initial radii. 

In all sets of models, the initial mass segregation allows 
to keep the majority of massive AGB stars bound in clusters 
with mass above 2 x 1O 3 M 0 . Without IMS, the MRR models 
for clusters with M c \ > 2 x 10 4 M o are more affected by the 
tidal field than other models because of their largest initial 
sizes, but the number of ejected stars remains fairly low in all 
cases. The fraction of massive AGB stars evolving in clusters 
is 68% (67%) with (without) IMS for the TRR models, and 
72% (75%) for the CRR models, respectively. These numbers 
are not very different from the value of 70% for the MRR 
models. 

The TRR models have a shallower slope of the mass- 
radius relation compared to the MRR models, 1/3 and 1/2, 
respectively. Therefore, the TRR models are characterised by 
lower intensity of intracluster UV radiation field in SCs with 
M c i < 2 x 10 4 M 0 , and more intense radiation field in clus¬ 
ters of higher mass. The former decrease in the field strength 
is more important, because about 80 % of progenitors of mas¬ 
sive AGB stars originate from clusters with M c \ < 2 x 10 4 M 0 . 
These results are quantified in Table 4 showing the fraction of 
massive AGB stars in which the UV field penetrates in the 
dust formation zone normalised to the number of these stars 
bound in clusters and to their total (bound+unbound) number. 
The numbers of AGB stars with the radii of SiO photodisso¬ 
ciation within the adopted values of R am d for the TRR mod¬ 
els are lower compared to the corresponding numbers for the 
MRR models (Tables 2 and 3), but the difference is at most 
10 %. 

The initial difference in the strength of the intracluster UV 
field compared the reference models is much stronger for the 
CRR case, although it is smoothed to some extent by the tidal 
disruption of low-mass clusters during the simulations. In 
contrast to low-mass clusters, the sizes of massive clusters 
increase by less than 2 times at most and preserve the initial 
difference, which is about 5 times between the CRR and MRR 


models for the 1O 5 M 0 star cluster. Impact of these compact 
massive SCs on the total fraction of stars affected by the UV 
field is diminished by the cluster mass distribution function, 
proportional to M d 2 (Sect. 2.2). Nevertheless, the photodisso¬ 
ciation radii of SiO are generally smaller in the CRR models, 
resulting in higher numbers of stars with Rpd.sio < Rmmi (Ta¬ 
ble 5). In particular, the fraction of bound massive AGB stars 
with Rpd.sio < 5 R* and D = 1CT 4 is 60% (42%) with(-out) 
IMS for the CRR models, which is more than two times higher 
than for the MRR models. The fraction of bound stars with 
SiO envelopes within 2 R t , our most conservative value of 
the condensation radius, is by 17% (25%) higher in the CRR 
models with(-out) IMS and D = 0 compared to the MRR. 

This small dependence on the adopted initial mass-radius 
relation is attributed to the important role of the galactic tidal 
field in evolution of low- and intermediate-mass clusters. It 
implies that our conclusions do not strongly depend on the 
initial conditions. 

5. DISCUSSION AND CONCLUDING REMARKS 

We investigate what fraction of massive AGB stars evolves 
in their parent star clusters and impact of cluster UV radia¬ 
tion on then' circumstellar environment. The field stars in our 
model originate (i) from low-mass clusters dissolve on the 
stage of emergence from the natal molecular clouds and (ii) 
by ejection from their natal star clusters. The latter process 
is modelled using N-body simulations of dynamical evolution 
of star clusters in the local Galaxy. We find that 70% of mas¬ 
sive AGB stars are members of their parent star clusters. This 
value can be lower if the fraction of stellar mass formed in star 
clusters is lower as suggested by some works (Piskunov et al. 
2008a; Kruijssen 2012). Half of all bound massive AGB stars 
resides in clusters with initial mass exceeding 5 x 10 3 M 0 . 

Intracluster UV radiation is 10 2 - 10 3 times stronger than 
the interstellar radiation field for stars 6-7 M 0 and gradu¬ 
ally decreases to a few times ISRF strength for 4 M 0 stars 
in all considered star clusters. The UV field strength shows 
10 times larger scatter towards lower values in clusters with¬ 
out initial mass segregation owing to more spatially extended 
stellar distribution. For this reason, the cluster UV photons on 
average penetrate deeper in the circumstellar shells of AGB 
stars in model clusters with initial mass segregation. 
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We find that on the initial stages of dust formation the clus¬ 
ter UV photons are able to dissociate the SiO molecules in 
the dust formation zone in a large fraction of massive AGB 
stars. This implies that dust condensation process may be 
suppressed or delayed by intracluster UV field, similarly to 
the effect of UV field of chromospheric origin (Beck et al. 
1992). The fraction of affected stars depends on the adopted 
initial mass-radius relation for model star clusters and loca¬ 
tion of dust-forming zone. The models with initial segregation 
of stellar mass are characterised by stronger UV fields and 
therefore smaller SiO envelopes. For our reference models, 
80% (60%) of AGB stars bound in star clusters with (with¬ 
out) IMS have the SiO photodissociation radii within 5 photo- 
spheric radii. Relatively to all (bound and field) massive AGB 
stars, this fraction constitutes 56% (41%). This value is even 
higher, 60% (56%), if we assume the constant initial radius of 
clusters supported by some observations (Larsen 2004). In¬ 
frared surveys for dust-enshrouded giants in young star clus¬ 
ters in Magellanic Clouds point to lack of reddened IR ob¬ 
jects (Van Loon et al. 2005b). However, they do find massive 
oxygen-rich AGB stars with dust emission signatures in some 
clusters. Our results imply that the sizes of SiO envelopes of 
AGB stars in cluster with ages < 200 Myr are smaller than 
those of their counterparts in the field population as seen in 
Fig. 9. This may be used as an indicator that an AGB star is a 
cluster member. 

Our findings have major implications for dust input from 
AGB stars at the solar galactocentric radius. Adopting the up¬ 
per limit for the extent of condensation zone R CO nd of 10 R, 
and assuming that dust formation is completely suppressed in 
stars with the SiO photodissociation radius R p d,sio < R C ond> 
we find that the effect of cluster environment on dust con¬ 
densation in CSE can reduce the total injection rates of sili¬ 
cate dust from AGB stars from 2.2 x 10 4 M Q pc 2 Gyr 1 to 
1.8 x 10 4 AL; pc 2 Gyr ', or at most by 20%. This value de¬ 
pends on the adopted dust yields and on the number ratio be¬ 
tween AGB stars of low and intermediate mass, which is de¬ 
termined by the star formation history. Applying the same as¬ 
sumptions to the model of stardust lifecycle in the solar neigh¬ 
borhood, we derived the lower limit for the revised contribu¬ 
tion of massive AGB stars to the presolar grain population 
of 40%. This value is 1.5 times lower than our earlier esti¬ 
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mate (Gail et al. 2009), but it is still too high to help explain 
non-detection of grains with isotopic signatures of hot bottom 
burning in meteorites. 

The lower mass limit for hot bottom burning adopted in the 
present work is 4 M 0 . If stars in mass range ~ 3 - 4 M 0 also 
experience hot bottom burning (Busso et al. 1999; Marigo 
2007; Ventura et al. 2012), they will be less affected by intr¬ 
acluster UV field, since after 220 Myr its strength becomes 
comparable to or weaker than the ISRF. Consequently, the 
fractions of massive AGB stars affected by the intracluster 
field would decrease, if we extend the considered mass range 
towards lower masses. 

In our simple estimate, we did not consider dissocia¬ 
tion of molecules which may be precursors of minor dust 
species with higher condensation temperature, e.g. corundum 
(ALO 3 ) and spinel (MgALO,}) grains. We demonstrated that 
a low dust-to-gas ratio of 10 4 may mitigate the penetration 
of the UV photons in the base of dust formation zone. Minor 
oxygen-bearing dust species with higher condensation tem¬ 
perature can thus assist formation of silicate grains by block¬ 
ing the UV irradiation. Clumpiness in the shell is another fac¬ 
tor presently neglected in our study which can provide local 
shielding from UV photons and enable dust formation. 

Most mass is lost by evolved stars during a few last ther¬ 
mal pulses of AGB evolution (superwind phase). A possible 
implication of the results derived in this work is that the su¬ 
perwind stage may be delayed in most massive AGB stars in 
clusters owing to suppression of initial condensation by UV 
irradiation of CSE. 
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